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Gate induced enhancement of spin-orbit coupling in dilute fluorinated graphene 
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We analyze the origin of spin-orbit coupling (SOC) in fluorinated graphene using Density Functional Theory 
(DFT) and a tight-binding model for the relevant orbitals. As it turns out, the dominant source of SOC is the 
atomic spin-orbit of fluorine adatoms and not the impurity induced SOC based on the distortion of the graphene 
plane as in hydrogenated graphene. More interestingly, our DFT calculations show that SOC is strongly affected 
by both the type and concentrations of the graphene’s carriers, being enhanced by electron doping and reduced 
by hole doping. This effect is due to the charge transfer to the fluorine adatom and the consequent change in 
the fluorine-carbon bonding. Our simple tight-binding model, that includes the SOC of the 2 p orbitals of F 
and effective parameters based on maximally localized Wannier functions, is able to account for the effect. The 
strong enhancement of the SOC induced by graphene doping opens the possibility to tune the spin relaxation in 
this material. 

PACS numbers: 73.22.Pr, 72.80.Vp, 85.75.-d,81.05.ue,71.70.Ej 


I. INTRODUCTION 

Spin based transport phenomena in graphene is a flourish¬ 
ing area of research due to the expected long spin relaxation 
lengths resulting from the small spin-orbit coupling (SOC) of 
carbon atoms and the high carriers mobility. Due to these ex¬ 
ceptional properties, the great potential of graphene for the 
study of new fundamental phenomena and for applications in 
nanoelectronics are now extended to the study of spin depen¬ 
dent phenomena with possible applications on spintronics. 1-4 
In this way the two-dimensional material, with unusual elec¬ 
tronic properties, became a model system in which the big 
challenge is to control and manipulate both, the charge and the 
spin degrees of freedom. Concerning the spin properties, the 
actual value of the spin relaxation rate of graphene carriers is 
still a controversial issue: while the spin relaxation times were 
initially thought to be very long, 5-9 different experiments 10-15 
suggest that they are much shorter than the theoretical predic¬ 
tions. This apparent controversy has recently been challenged 
by the way experiments have been interpreted. 16 

Whatever the source of the spin relaxation mechanism, 
an important issue for spintronics applications is the abil¬ 
ity to control it. 3,4 Sources of spin relaxation in graphene 
are usually attributed to magnetic defects (vacancies 17-22 or 
adatoms 9,23-37 ), to the presence of extrinsic (adatom induced) 
Soc 9 ’ 38 ripples 39 or, more recently, to a combined effect of 
spin-orbit and pseudo-spin physics 40 

Adatoms, in particular, provide a possible course to engi¬ 
neer spin based effects. The nature of the SOC, however, 
might be different for different adatoms. Very light atoms, 
like H, which are adsorbed on top of a single C atom and have 
a small intrinsic atomic SOC, introduce SOC to the graphene 
carriers by distorting the otherwise flat graphene sheet. 9 This 
is due to the sp 3 like structure adopted by the hybridized C 
atom that induces a local coupling between the carbons p z 
and cr orbitals. In that case, the effect of the atomic SOC 
of the C atoms changes (locally) from being a second order 
effect to a first order one with the corresponding increase of 
the effective spin-flip processes in the graphene carriers. A 


characteristic of this mechanism is that the resulting SOC is 
proportional to the local lattice distortion. 9 On the other hand, 
more heavy atoms can introduce spin effects by processes that 
involve their own intrinsic SOC, i.e. where the spin-flip oc¬ 
curs in the adatoms orbitals and not on the C atoms. In the 
general case both mechanisms are present. 

The fluorine atom on graphene is adsorbed in a top position. 
It bounds covalently to the C atom below it (we will refer to 
it as the Co atom), therefore introducing a local distortion in 
the graphene lattice. However, owing to its strong electro¬ 
negativity, there is a charge transfer from the graphene sheet 
to the F adatom. 36,41,42 This charge transfer to the F adatom 
can be controlled by changing the carriers concentration of 
graphene. Associated to this, there is a sp 3 -sp 2 like crossover 
of the hybridization of the Co atom 41,42 Since this crossover 
causes a strong reduction of the lattice distortion, one would 
expect, as suggested by the scenario described in Ref. [9], that 
the induced SOC associated with it will also decrease, thus 
providing a controllable way to reduce the spin relaxation. 

Here we show, using Density Functional Theory (DFT) cal¬ 
culations, that, contrary to this expectation, the SOC in dilute 
fluorinated graphene is strongly enhanced by electron dop¬ 
ing. This is due to the fact that the main contribution to the 
SOC comes from the intrinsic SOC of the F adatom 43 and not 
from the C atoms as the lattice distortion mechanism requires. 
This result is validated by a simple tight-binding model, con¬ 
structed using the Wannier functions parameter obtained from 
the ab initio calculations, which accounts both qualitatively 
and quantitatively for this effect. 

The rest of the work is organized as follows: DFT re¬ 
sults for the band structure and the projected density of states 
(PDOS) for different geometries are presented in Sec. II. A 
tight-binding model based on the calculations of the maxi¬ 
mally localized Wannier is introduced in Sec. Ill to discuss 
the microscopic mechanism leading to SOC. We finally con¬ 
clude in Sec. IV. 
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FIG. 1. (a) Electronic band structure for a 7 x 7 supercell containing 98 C atoms and a single F adatom. The horizontal dashed line indicates 
the Fermi Energy C F which is set to zero, (b-e) projected density of states (PDOS) onto the p z (solid line) and p^ orbitals (dashed line) for 
the F adatom, the carbon bound to it (C 0 ), one of the three nearest C atoms (C„) of the latter, and a ‘bulk’ (far from the adatom ) carbon atom 
(C b ). a , p and y label the closest bands to C F . 



II. DFT RESULTS 

The DFT calculations were performed with the Quantum 
Espresso package 44 employing density functional theory and 
the Perdew-Burke-Ernzerhoff (PBE) exchange-correlation 
functional 45 An ultra soft description of the ion-electron in¬ 
teraction was used 46 together with a plane-wave basis set for 
the electronic wave functions and the charge density, with en¬ 
ergy cutoffs of 70 and 420 Ry, respectively. The electronic 
Brillouin zone integration was sampled with an uniform k- 
point mesh (3x3x1) and a Gaussian smearing of 0.005 Ry. 
The two-dimensional behavior of graphene was simulated by 
adding a vacuum region of 12 A above it. All the structures 
were relaxed using a criteria of forces and stresses on atoms 
of 0.005eV/A and 0.3GPa, respectively. The convergence 
tolerance of energy was set to 10 -5 Ha (1 Ha = 27.21 eV). 
To correct for the dipole moment generated in the cell and 
to improve convergence with respect to the periodic cell size, 
monopole and dipole corrections were considered 47 This is 
particularly important in the doped cases. Doping of the unit 
cell (added/removed electrons) where compensated by an uni¬ 
formly distributed background charge. 

A. Electronic band structure 

Let us first analyze the band structure of fluorinated 
graphene in the absence of additional doping (neutral case). 
Since we are interested in describing the SOC introduced 
by a single impurity (corresponding to a diluted fluorinated 
graphene sample), for our DFT calculations we use super cells 
as large as posible (limited by the computational resources). 
The use of large super cells is also important in the case of 
F to avoid long range Coulomb interactions among impu¬ 
rities as F adatoms acquire some charge when absorbed on 
graphene 41-43 


Figure 1(a) shows the DFT electronic band structure for a 
7x7 supercell containing 98 C atoms and a single F atom. 
The Fermi energy E p has been set to zero, E? = 0. The path in 
the reciprocal space is labeled using the standard notation for 
the hexagonal Brillouin zone of the supercell lattice. Though 
negligible on this scale, each band is split in two due to the 
SOC—this is analyzed in detail in the next section. Since one 
would naively expect to observe graphene’s Dirac cone at the 
K point, it is instructive to mention why that is not the case: 
our DFT supercell contains a single F atom and consequently 
all impurities in our periodic system are on the same graphene 
sublattice. This breaks off the sublattice symmetry and opens 
a gap in the graphene spectrum. 48,49 Such gap is proportional 
to the impurity concentration and hence to the supercell size. 
While this effect should not have any relevance for the cal¬ 
culation of local parameters for a large enough supercell, it 
is important to keep it in mind for a proper interpretation of 
the results—the breaking of the sublattice symmetry can be 
eliminated from the calculations, at the price of increasing the 
F concentration or increasing the supercell size, by including 
two F atoms on the supercell, one on each sublattice. 

To analyze the character of the bands, it is useful to look at 
the weight of the corresponding states on the atomic orbitals 
of each type of atom. That is, to look at the projected den¬ 
sity of states (PDOS). This is done in Fig. 1 for four different 
atoms in the supercell: F (b), Co (c), one of the three C n car¬ 
bon atoms around Co (d), and a ‘bulk’ (far from the adatom) 
C atom (e). The solid line corresponds to the p z orbitals while 
the dashed one corresponds to the p x (y) orbitals—in this en¬ 
ergy range the weight on the s orbitals is smaller than 0.02 
eV -1 . It is apparent from the figure that the central band (la¬ 
belled by P) and the band below it (y) are the ones that con¬ 
tain the impurity states as they carry a large weight on the p z 
orbitals of the F adatom, the Co atom, and its three nearest 
neighbors C n atoms. 37 

With the aim of understanding the adatom induced SOC, a 
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FIG. 2. Same as Fig. 1 but for a 6 x 6 supercell containing 72 C atoms an a single F adatom. The last rightmost panels correspond to a C atom 
far from the adatom in the same (e) and the opposite (f) sublattice. Notice that in this case there is a clear Dirac-like band and that the F and 
the C 0 atoms are essentially decoupled from it. 


very interesting situation occurs for supercells with a single F 
and sizes of the form 3 n x 3 n with n an integer number. In this 
case, the perturbation potential induced by the adatom U(r) 9 
assumed to be local, has the right symmetry as to mix the two 
non-equivalent graphene’s Dirac cones. Namely, the Fourier 
transform of the periodic adatom potential contains non-zero 
components U & with k connecting the two cones. Such ad¬ 
mixture leads to a partial decoupling of the graphene n band 
from the adatoms. This closes the gap induced by the sub¬ 
lattice symmetry break and restores a single Dirac cone. In 
addition, the impurity band is more clearly developed. The 
corresponding band structure is presented in Fig. 2 for a 6 x 6 
supercell, together with the PDOS. Notice that a Dirac cone 
can be clearly identified and that in those bands there is essen¬ 
tially no weight of the fluorine’s p z orbital which, in this en¬ 
ergy range, is mainly concentrated on the rather flat impurity 
band (the band in Fig. 2(a)). This decoupling, that occurs for 
a particular geometry, and though rather artificial, will help us 
to separate different contributions to the SO splitting of the 
bands. 


B. Gate dependence of the spin-orbit splitting. 

We now turn our attention to the splitting of the bands in¬ 
duced by the SOC and, in particular, to the effect of doping on 
such splitting. 

The main contribution to the SO splitting in this system 
arises from the atomic SOC of the C and the F atoms. The 
contribution from the C atoms is known to be a second or¬ 
der effect in flat graphene own to the reflexion symmetry of 
the graphene’s plane: symmetry prohibit a direct coupling be¬ 
tween the p z and p*^ orbitals of adjacent C atoms—recall that 
atomic SOC mixes p orbitals on the same lattice site. This 
symmetry is locally broken in the presence of an adatom that 
sits on top a C atom due to the lattice distortion it introduces. 
This was shown 9 to enhance the SOC by breaking the above 
mentioned selection rule. The resulting SOC is proportional 
to the distortion of the lattice. The latter depends, in the case 
of F adatoms, on the doping level of the graphene sheet due 
to the charge transfer from the graphene to the F adatom. 37,42 


Hence, one expects a strong dependence of the C contribution 
to the SO splitting upon doping. Namely, since doping with 
electrons (holes) reduces (enhances) the sp 3 character of the 
local hybridization of the Co atom, we expect this contribution 
to the SO splitting to decrease (increase). 

The other source of SOC is the adatom itself. Contrary to 
the case of H, the SOC in the F atom is not negligible and must 
be accounted for. How this contribution changes with doping, 
if it does, is however not obvious a priori. 

The DFT results for the SO splitting of three bands for the 
case of the 7 x7 supercell (indicated as a , /?, and y in Fig. 1(a)) 
is shown in Figs. 3(a-c) for five different doping concentra¬ 
tions: Sn = 0, ±1/2, ±1 is the number of additional electrons 
per unit cell (hence Sn = 0 corresponds to the neutral case). 
We find that the SO splitting changes for all the bands, be¬ 
ing enhanced by electron doping and reduced by hole doping. 
This is just the opposite behavior that is expected based on the 
distortion induced SOC. We must then conclude that the main 
contribution to the observed splitting is not the atomic SO of 
the C atoms but the one coming from the F adatoms. 43 To 
verify this, we repeated our calculation but without including 
the SOC on the F—this is done by using a scalar relativistic 
pseudo-potential for the F atom as implemented in the Quan¬ 
tum Espresso code 44 ). The results are shown in Figs. 3(d-f). 
There are two important points to notice: (i) the magnitude 
of the splitting is reduced by a factor of approximately five, 
consistent with the difference in magnitude of the atomic SO 
between F and C (~ 50 and ~ 10 meV, respectively), and sig¬ 
naling that the main source of SO is absent; (ii) the SO split¬ 
ting shows now the expected behavior as a function of doping 
for a deformation induced SO: it is reduced by electron doping 
and enhanced by hole doping. 

To further verify the origin of the SO splitting in fluorinated 
graphene we have done the calculation keeping the SOC in the 
F atom but using a 6 x 6 supercell that leads to the decoupling 
of the 7i band from the adatoms as discussed in the previous 
section. Quite interestingly, the bands show now different be¬ 
haviors (Figs. 4(a-c)): those with a small weight on the flu¬ 
orine’s p z orbitals (bands a and y) have a SO splitting com¬ 
patible in magnitude with a deformation induced SO and the 
corresponding dependence with doping. On the other hand, 
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FIG. 3. (color online) SO splitting of the a, ft, and y bands indicated 
in Fig. 1(a) for the cases where the SOC of the F atom is included in 
the DFT calculation [(a-c)], and those where it is removed [(d-f)] as a 
function of the doping, 6n = 0, +1/2, +1 electrons per unit cell. The 
arrows indicate the direction of increment of 6n while the multipli¬ 
cation factors represent the scaling used in the energy scale. Notice 
that the doping dependence is just the oposite with and without SOC 
in the F atom. 


the band with a large weight on the F (J3 ) shows a large SO 
splitting and the opposite doping dependence. Again, this is 
confirmed by removing the SOC in F as shown in Figs. 4(d-f). 


III. TIGHT-BINDING MODEL 

The above results clearly demonstrate that in fluorinated 
graphene the main source of SOC are the F adatoms, in agree¬ 
ment with Ref. [43]. However, the origin of its strong de¬ 
pendence on doping is not at all clear. To better understand 
the underlying microscopic mechanism we now construct a 
single particle tight-binding model that includes the most rel¬ 
evant orbitals, and later on the atomic SOC of F. Namely, we 
consider a small cluster consisting of one F and four C atoms 
(see Fig. 5) embedded in graphene. The Hamiltonian is given 
by *H = + ‘Hfc- The first term describes the isolated 

F atom 


= Yj S P P lo-Pt;o- + Y £ ' ,S 'A > (U 

£,cr cr 

where creates an electron with spin cr in the 2p^ orbital 

(£ = x, y, z) and vj- creates an electron with spin cr in the 2s 
orbital. Given the symmetry of the system, it is convenient, 


FIG. 4. (color online) SO splitting as in Fig. 3 but for the 6 x 6 su¬ 
percell. Notice that only the impurity band (/?) maintains the doping 
dependence of the F induced SO splitting, the rest behave as the C 
induced SOC near the Y point. 


for the definition of the hopping parameters and for the sake 
of comparison with the Wannier functions described below, to 
work in a hybridized basis. To this end we define the following 
creation operators, 9 


fl=Asl+^A^pl, 

Vl-A 2 t A t 

, t _ t A t 

fl ° V 3 ,<r 

Vl-A 2 t A t 


+ 



1 

l 


l 

vt 

1 


Pycr ’ 




( 2 ) 


The parameter A serves to interpolate between two extreme 
cases: (i) A = 0, in which case = p\ a while the other or¬ 
bitals, correspond to the standard sp 1 hybrid orbitals; (ii) 
A = 1/2, where all orbitals correspond to the sp 3 hybridiza¬ 
tion, with one of them pointing in the z direction. In terms of 
these orbitals, the fluorine Hamiltonian takes the form 


= Y Ez + Y E f flrfbr + Y h f Afjcr 

<j i,cr i± j,CT 

+ X t2f {fir fUr + fkrfzo) ’ 
i,cr 


( 3 ) 
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FIG. 5. (color online) Schematic representation of the hybridized 
orbitals used in the tight-binding model with their hopping matrix 
elements. 


where i,j- 1,2,3, and 


. x (l-A 2 ) ( 2 + A2 ) 

£ -~ = ( 1 “ a2 ) e p + a2£s > e f = 1 £s + 1 E P ’ 




( 1-a2 )/„ . _ aVi-a 2 ^ ^ 

^ ( e i E p) . hf - _ ( e i E p) ■ (4) 


The Hamiltonian of the C atoms is 


Kc = J]eo cl kr C { , T + J]e l c] a c ilr + /> + c l c 0o-) 


+t 2 Z 4 T C jcr + 'Hg , 
i±j,<r 


(5) 


here c and (/ = 1,2, 3) create electrons with spin cr at 
the p z orbitals of the central C atom (Co) and the side carbon 
atoms (C„), respectively, and Ei g is the Hamiltonian of the rest 
of the p z orbitals of the graphene sheet with energy s g and 
hopping t. Finally, TTpc includes the hybridization between 
the F and C orbitals, 

■Hfc = V Z [f^c 0(r + clj zo ) + V'Yj {fzo- c io- + c lfzp) (O 
cr icr 

+ W Z {fl c 0 <r + 4 jbr) + Z W ‘j (& C *r + ’ 

i,a ij,cr 


where takes only two values: w\ if the f? orbital ‘points’ 
towards the ct orbital and W 2 otherwise (see Fig. 5). 

The next step is to properly estimate all the parameters: en¬ 
ergies of the hybridized orbitals and hopping matrix elements. 
This can be achieved with the help of the Wannier90 code 50 
as shown in the next section. 


A. Wannier functions 


Doping 

-1 

1 

2 

neutral 

+ 2 

+ 1 

£ z 

-10.43 

-9.63 

-8.72 

-7.50 

-5.40 

£ f 

-12.29 

-11.45 

-10.35 

-8.86 

-6.27 

£0 

-0.08 

0.25 

0.38 

0.31 

0.36 

Si 

-2.55 

-1.91 

-1.30 

-0.80 

-0.27 

f 

-2.45 

-2.37 

-2.31 

-2.34 

-2.71 

*2 

-0.07 

-0.09 

-0.11 

-0.09 

0.13 

b/ 

-5.34 

-5.25 

-5.15 

-5.02 

-4.85 

t 2 f 

-3.99 

-4.07 

-4.17 

-4.25 

-4.33 

V 

3.34 

3.31 

3.11 

2.48 

1.11 

V' 

-0.21 

-0.18 

-0.14 

-0.06 

0.09 

w 

-3.47 

-3.43 

-3.18 

-2.42 

-0.88 

W\ 

-0.54 

-0.56 

-0.56 

-0.53 

-0.42 

w 2 

0.30 

0.28 

0.25 

0.20 

0.07 

s g 

-2.08 

-1.53 

-0.99 

-0.52 

-0.34 

t 

-2.85 

-2.86 

-2.87 

-2.88 

-2.90 


TABLE I. Tight-binding parameters (in eV) calculated by the Wan- 
nier90 program for different doping concentrations (number of addi¬ 
tional electrons per unit cell) in a 4 x 4 supercell. For a description 
of the parameters see Fig. 5. 


only approximately, depending on the number of orbitals used 
in the calculations. To this end, we use the Wannier90 code 
in a 4 x 4 supercell. While our calculations include the cr 
orbitals between the C atoms, required to properly describe 
the DFT band structure, here we present only those parameters 
that are relevant for the simplified tigh-binding Hamiltonian 
described above which is enough to capture the SO splitting of 
the energy bands around E F . Table I presents such parameters 
for several doping configurations. 

The maximally localized Wannier functions obtained for 
the Cq, cj, f z , and orbitals are shown in Fig. 6. It is ap¬ 
parent from the figure that the F orbitals can be interpreted as 
the hybridized atomic orbitals described above while the Co 
and C n orbitals need to be considered as ‘effective’ orbitals 
(the former being more sp 3 -like and the latter as a tilted p z - 
like one). 

The interpretation of the F orbitals as hybridized atomic or¬ 
bitals is also supported by the following: our tight-binding 
model has three parameters (s S9 s p , A) from where four pa- 


Doping 

e, [eV] 

S P [eV] 

A 

qf [eV] 

/.calc I + 

t 2 f E 2 f 

-1 

-26.45 

-6.95 

0.42 

-4.31 

1.08 

1 

2 

-25.38 

-6.20 

0.42 

-4.24 

1.04 

neutral 

-24.17 

-5.20 

0.43 

-4.26 

1.02 

+ ! 

-22.56 

-3.84 

0.44 

-4.29 

1.01 

+ 1 

-19.95 

-1.42 

0.46 

-4.39 

1.01 


A simple way to build a tight-binding model from the DFT 
results is to find the maximally localized Wannier functions 51 
that describe the DFT band exactly on the energy range of 
interest (around E F , for instance) and the rest of the spectrum 


TABLE II. Parameters of the model Hamiltonian s s , s p , and 
A are calculated by fitting the values of s z , Sf, and t\f obtained by 
the Wannier functions method. The value of ty C is derived from the 
former, being consistent with the Wannier calculation. 
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FIG. 6. (color online) Maximally localized Wannier functions as ob¬ 
tained with the Wannier90 code and used for the tight-binding model. 
They correspond to: (a) cj, (b) one of the cj orbitals, (c) fj, and (d) 
one of the f? orbitals. 


rameters are obtained [cf. Eqs. (4)]. These four parameters 
are calculated independently in the Wannier90. Hence, we can 
use three of them, say (s z , Sf, tif), to calculate ( e s , s p . A) and 
from them the value of ty C to be compared with the value ty 
shown in Table I. Such comparison, together with the values 
obtained for s s , s p , and A is shown in Table II. Our results 
show that, within this model, the main effect of doping on the 
F’s orbitals is to shift their energy up as more electrons are 
added to the supercell and to slightly change the hybridization 
parameter A (< 10%). This energy shifting is to be expected 
from a naive mean field argument as doping induces a charg¬ 
ing of the F adatom. 41,42 We would like to point out that while 
with this very simple model one cannot expect to fit the value 
of t 2 f, the calculated hopping term (ty C ) is quite reasonable 
and follow the overall trend with doping—except for the ex¬ 
treme hole doping case. 

It is worth noticing that while there is an overall shift of all 
energies with doping, the important magnitude is their relative 
shift, which is by far larger in the F atom. This is the reason 
for the change of the SOC as we discuss below. 

Finally, we emphasize that the present calculation was done 
in a 4 x 4 supercell which is relatively small as to avoid the 
finite size effects related to the long range Coulomb repulsion 
between F adatoms 42,43 —in our case the size is limited by the 
heavy computational resources needed for the calculation of 
the Wannier orbitals. Therefore, we may expect that our tight- 
binding parameters will not be reliable for large electron dop¬ 
ing, where the finite size effect are stronger. Yet, we will show 
that they provide a good description of the SO splitting of the 
bands. 



FIG. 7. (color online) Comparison of the SO splitting of the a,j3 , and 
y bands obtained with the simplified tight-binding model (solid lines) 
and the DFT results of Fig. 3 (dashed lines), excluding the case of 
+1 electron doping. To obtain a good agreement we use a = 40meV 
while the rest of the tight-binding parameters are taken from Table I. 


B. Spin-orbit coupling 


As our DFT results show that the main source of SOC 
comes from the F adatom, we will include in our tight-binding 
model only the atomic SOC of the F—an extension to include 
the effect of the SOC on the C atoms is straightforward. Such 
term is of the form p° = a L • S, with L and S the orbital 
angular momentum and spin operators of electrons in the p 
orbitals, respectively, and the SO parameter a ~ 50meV for 
F. 43 ’ 52 In the {p* T , p vT , p- T , p xl , p vi , p 4 } basis it takes the fol¬ 
lowing matrix form 53,54 


r 0 -i 0 0 0 1 ' 

i 0 0 0 0 -/ 

<j_^so _ ® 00 0 - 1 / 0 

F “ 2 0 0 -1 0 /0 
0 0 —i -i 0 0 
A i 0 000, 


( 7 ) 
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FIG. 8. (color online) Scheme showing the most relevant hopping 
parameters of Hamiltonian 77. The orbital states decoupled by sym¬ 
metry only couple through the spin flip-term of the fluorine’s SOC. 


By adding 77p° to the tight-binding Hamiltonian presented in 
the previous section we can now compute the splitting of the 
a , ft and y bands for a given supercell and compare it with 
the DFT data. This is done in Fig. 7 using the parameters of 
Table I and taking a ~ 40 - 50meV. There are several points 
to emphasize: (i) there is an overall good qualitative agree¬ 
ment, which further corroborates our model; (ii) the value of 
the splitting and its behavior with doping is the correct one, 
indicating that the F adatom is the dominant source of SOC; 
(iii) our model overestimate the effect of electron doping (the 
+1 case it is strongly overestimated hence is not shown in the 
figure). Presumably this is due to the finite size effect dis¬ 
cussed in the previous section that affects our tight-binding 
parameters; (iv) the central band (J3 ) is the one that is better 
described, as could be expected. 

To better understand what are the microscopic spin-flip pro¬ 
cesses that lead to the splitting of the bands it is useful to 
take full advantage of the symmetry of the system around the 
adatom and define new symmetrized orbitals as follows, 


f »• = _ y ft 

-'Acr J icr ’ 


'• + = — Yd, (8) 


V3 tr 


/ j icr ’ 


and 

f t _ PAzAl11Al 
/+cr V3 


P C lr ~ C L + P* C 


t _ H Icr ^2cr 


3 cr 


V3 


/_cr V3 


«2 t , t ,o*2 t 
t _ Z 3 L l<r _ c 2er L 3cr 

V3 


( 9 ) 


where ft - exp(/7r/3). Note that fl a = Vl - A 2 si - A p\ a 
is the combination of the s and p z orbitals that is orthogonal 
to f z(r . In this basis, the full Hamiltonian decouples into or¬ 
thogonal subspaces with different spatial symmetry. The key 
point is that this orthogonal spaces only couple through the 
splin-flip term of the fluorine’s SOC. The Hamiltonian of the 
rest of the graphene sheet can also be separated in two sectors 
with different spatial symmetry—each one of these couple to 
one of the c\ or c\ orbitals. 

A scheme of the hierarchy of the hoppings matrix elements 
in 77 is shown in Fig. 8 for one of the two spin sectors— 
the other one is analogous but with the time reversal part¬ 
ners. Since first order (direct) spin-flip processes only occur 
between the f zcr , fl CT and the f ±& orbitals of the F adatom, 
any effective SOC between C atoms must be induced by vir¬ 
tual processes that involve them at intermediate steps. This is 
the way a spin-splitting is induced in all bands. The scheme 
makes evident a couple of things: (i) no virtual spin-flip pro¬ 
cess is allowed between the F and the Co atoms as they belong 
to the same spin subspace, in agreement with the analysis of 
Ref. [43]; (ii) virtual spin-flip processes between C atoms, 
necessarily involve the mixing of the subspaces with oppo¬ 
site spatial symmetry and hence the coupling between the fl a 
and the c\ a orbitals is crucial. The effect of doping enters by 
changing (reducing or increasing) the energy difference be¬ 
tween the F and the C orbitals, making the virtual spin-flip 
processes more or less effective. 

Finally, it is useful to estimate from our model the effective 
spin-flip coupling between the C n atoms around Co (ie, an ef¬ 
fective term of the form A c\^c 2 ^ + A* for instance) in 

order to compare with other models in the literature. 9,43 Fol¬ 
lowing the standard procedure to eliminate the intermediate 
states using Green functions, we get A ~ 2meV for the neutral 
case. This effective SOC changes with doping: it decreases up 
to roughly 20% for hole doping and increases up to 60% for 
the cases shown in Fig. 7. This parameter is related to the 
parameter Ap IA introduced in Ref. [43]. Our result has the 
same order of magnitude, indicating a large SOC induced by 
the adatom, but differs in a factor ~ 3. 


IV. SUMMARY 


We have analyzed the mechanisms governing the appear¬ 
ance of SOC in fluorinated graghene. Our results allow to 
disentangle the different contributions of the SOC induced on 
the graphene carriers and to make predictions on its evolu¬ 
tion with electron and hole doping. DFT calculations of a F 
atom adsorbed on different graphene supercells show signif¬ 
icant spin splittings of the low energy bands and give a first 
indication of the dominance of the atomic SO of the adatom 
in this case. The DFT results, combined with the Wannier90 
code, allow finding the maximally localized Wannier func¬ 
tions that correctly describe the bands on the energy range 
of interest to build an effective tight-binding model that in¬ 
cludes effective 7r-orbitals of the C atoms and the 2s and 2 p 
orbitals of the F atom. The so obtained tight-binding model 
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fits the spin-splitting of the low energy bands and allow for a 
microscopic interpretation of the origin of the SO effects and 
its doping dependence. We conclude that: The SOC in fluo- 
rinated graphene is dominated by the intrinsic SO interaction 
in the F atoms, in agreement with the results of Ref [43]. The 
charge transfer from the graphene to the F orbitals has an im¬ 
portant effect on the final spin-flip coupling of the graphene 
carriers. This is mainly due to the charge transfer induced en¬ 
ergy shifts of the fluorine sp orbitals that are mixed by the 
SOC (orbitals with energies e z and 6/ in Hamiltonian (3)). As 
this charge transfer can be controlled by doping, the final SO 
effects can be controlled by gating the sample. While hole 
doping produces a small decrease of the SOC, electron dop¬ 
ing can produce a significant increase of the effect. 

It is worth noting that our results show that the spin-splitting 
at the K point of the Brillouin zone is always extremely small. 
This is due to the symmetries of the supercells used in the 
calculation and by no means implies that SO effects on the 
low energy states of samples with a random distribution of F 
impurities will be negligible. To see this, consider the single 
impurity case of the system described by Hamiltonian of Eq. 
(??). There, the important SO parameter is the atomic SOC 


of the F atom a. As discussed in the previous section, this pa¬ 
rameter leads to an effective spin-flip coupling between the C 
atoms, an effect that can be described by an effective Hamil¬ 
tonian including spin-flip processes in the graghene carriers 
around the impurity site. 43 The spin-flip cross section for this 
effective impurity model was calculated in Ref. [9]. An im¬ 
portant result of our work is to show that, for fluorine impuri¬ 
ties, the spin-flip coupling between the C n atoms ( A) is much 
larger than what we could expect from the atomic SOC of the 
C atoms and, even more important, that its magnitude can be 
controlled by gating the sample. 
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